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Part  I 

FORWARD 

A carefully  developed  body  of  knowledge  concerning  estimation  in 
sonar  detection  and  classification  has  been  developed  and  made  available  to 
Naval  Ships  Syn terns  Command  in  previous  reports.  In  section  I ve  summarize 
definitions  of  measurement  space,  feature  space,  decision  space,  and  the 
decision  rule,  and  we  comment  on  the  generalised  k nearest  neighbor  decision 
rule,  nonparametric  feature  selection,  cluster  mapping,  and  the  Rayes  solution. 
Then,  in  Chapter  II,  we  introduce  the  new  appronch  to  pattern  recognition 
which  provides  for  inserting  problem  knowledge  into  the  decision  making  pro- 
cess. This  approach,  although  not  so  new  in  "idea"  is  very  new  in  "action". 

We  discuss  nonlinear  mappings  we  have  discovered  which  can  be  used  to  insert 
a priori  problem  knowledge.  We  show  how  features  can  be  designed  first  fen 
submarine  targets  and  then  for  specific  kinds  of  nonsubmarine  targets,  lly 
using  a priori  problem  knowledge,  pattern  recogniti on  can  t nko  pi aeo  with 
no  training  samples;  the  training  samples  then  add  ice  cream  to  the  pie. 

Section  IV  provides  an  introduction  to  performance  results  after  pro- 
cessing Rogers  data. 

I Theoretical  Foundation 
A.  Introduct ion 

Appropriate  definitions,  notation,  and  operations  have  been  fairly  well 
established  for  computerized  detection  and  classification.  In  this  section  we 
7rovlde  a tutorial  and  indicate  where  improvements  can  be  jt  do. 

Measurement  or  Observation  Space 

Measurements  of  n sonar  echo  or  sequence  of  sonar  echoes  provide  components 
of  a vector  jk,  “ , x2»,,,»xiJ  called  the  observation  vector.  These  meajuremonts 

•re  made  on  the  audio  waveform,  video  scan  waveform,  along  with  associated 


i 


b. 


L 


measurements  of  target  range  R,  true  tearing  0,  and  target  track  aspect  ljf. 
feature  Space 

The  process  of  obtaining  features  from  the  observation  vector involves 
inserting  problem  knowledge.  Features  important  for  detecting  a submarine  are 
different  from  features  important  for  detecting  a noncubmarine.  For  example,  a 
nonsubmarine  target  may  consist  of  multiple  targets  in  range,  bearing  scan,  or 
both;  features  which  measure  these  properties  are  thus  important  for  detecting 
nonsubmarines.  On  the  other  hand,  an  estimate’ target  length  is  an  important 
feature  for  detecting  submarine  targets.  It  is  important  to  provide  that  features 
useful  in  detecting  one  type  of  tarret  are  not  a nusicance  in  detectin'*  another  type 
of  target. 


Decision  Space 

Simply,  the  decision  space  consists  of  M classes  with  corresponding  a 
priori  probabilities  P^,  i - 1,2, ...,M.  We  may  consider  the  decision  space 
one  dimensional  with  K points. 

The  Uninn e Decision  Rule 

Loosely  speaking,  the  unique  decision  rule  is  a Bayes,  minimum  conditional 

risk  rule  where  the  density  f(x|i)  and  a priori  class  probability  P.  are  calculated 

for  all  M classes.  It  is  necessary  to  learn  f(x:|i)  using  both  Problem  knowledge 

and  training  vectors  x^,  x^>***>  x1  for  this  i*'*1  class.  The  density  f(x|i)  is 

X — n^  1 

characterized  by  a set  of  parameters  where  Jm  is  related  to  the  features  for 
the  class. 

£±_JLbe.  New  Prarrmatlc  Approach  to  Pattern  Recognition 

We  at  Purdue  University  have  made  several  contributions  to  pattern  recognition 
including: 

The  Generalized  K-Nearest  Neighbor  Decision  Rule,  published  in  Information  and 
Control 

Kon-Pa ram et ri c Feature  Selection,  published  in  IEEE  Transactions  on  Information 
Theory 


lAtop. 


f 
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Cluster  Mapping,  with  Experimental  Computer  Graphics,  published  in  IEEE 
transactions  on  Computers. 

Baves  Solutions  to  Pattern  Recognition,  published  in  IEEE  transactions 
on  Information  Theory. 

how,  however,  we  are  working  On  a more  advanced  approach  to  pattern  recognition 
which  provides  for  inserting  problem  knowledge  along  with  training  samples  into 
the  automation  process.  Subsequent  sections  will  describe  the  proceedure. 

II  The  Kevr  Proceedure 

The  approach  is  to  use  problem  knowledge  to  transform  observation  vectors 
jx*  from  class  i to  a lower  dimensional  space  where  the  features  are  uncorrelated. 

A nonlinear  transformation  must  be  constructed  for  each  class;  this  nonlinear 
transformation  may  be  thought  of  as  resulting  ultimately  in  a nonlinear 
matched  filter  vjhich  can  be  updated  using  training  samples.  The  proceedure  is 
best  described  in  terms  of  the  following  operations  used  to  construct  the  nonlinear 
transformations . 

Within  Subset 

This  operation  selects  measurements  x , x , x from  the  sth  sonar  ping 

S1  s2  SL 

to  which  nonlinear  operations  are  applied.  These  measurements  are  in  the  sth 

vector 

Between  Subset 

This  operation  selects  measurements  from  each  of  a sequence  of  v pings 
iSi*  7^2*  * • • which  nonlinear  operations  are  applied.  One  such  nonlinear 

operation  is  the  "Starlight"  type  operation. 

Pcp.ion 

This  operation  selections  a region  of  the  observation  space  V^.  Any  points 
in  this  region  are  processed  with  a sequence  of  operations  specialized  to  samples 
restricted  to  this  region. 
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Con-plat i op.  In^'j-tcr 

, .*  linea**  or  nonlinear  function  of  neap  -ements  selected  by  operations 
"vithin  subset"  or  "between  subset"  gives  rise  to  a new  feature.  The  following 
is  a list  of  nonlinear  functions  which  nay  or  my  not  be  directly  used  in 
sonar  problem: 

ratio 

sum 

weighted  sum 

blanking 

quadratic 

ring 

product 

subgroup 

More  specialized  nonlinear  relationships  have  been  found  from  sonar  problem 
knowledge  and  are  described  next. 

AJVjd  min?  ho.nl i near  Filters  for  Detecting  Subs 

Measurement  vectors^,  x^,...,x^  are  taken  for  v successive  pings  with 

AAA  A A A A 

Jp  B [ ®r»  ^r3»  r 1,2, ... ,v 

A A A 

where  t is  an  estimate  of  target  length,  E of  target  echo  energy,  HP  of 

A A A 

noise  energy,  ^ of  aspect  angle,  Rp  of  target  range,  C of  target  bearing, 

A 

and  Nr  of  number  of  target  highlights.  When  the  same  class  of  target  is  active 
for  all  v echos,  the  sequence  vector  of  measurements 

= x2,...,xy] 

is  constructed.  If  there  are  L entries  in  x then  there  are  bv  entries  in 

~r 

xy  and  they  can't  be  uncorrelatcd.  V’e  know  that  target  length  6 is  a fixed 
parameter  and  it  is  reasonable  to  use  the  sum  operation. 


1 'T 

£ " *■  ; <v  i 
v L,  r r 


where  <y  r'  E INI’  is  an  estimate  of  signal/noise  ratio  for  the  r^^  echo. 

r r*  r 

Th»  target  reflectively  is  a reasonable  feature  characterizing  submarine 
targets.  An  estimator  for  target  reflectivity  o is  signal/noise  ratio  normalized 


for  range. 


2 IK 

( £g*)  c')2  - f- 

Ih 

r*l 


whore 


v L.  r 
r~  1 


The  number  of  highlights  is  another  reasonable  feature  characterizing 
submarine  targets.  Since  the  number  of  highlights  in  a submarine  echo  deponds 
on  the  submarine  aspect  angle,  an  estimator  for  the  number  of  highlights  is 


» “ i £ "r  c°=  *, 


vhere  N is  a measurement  of  the  number  of  peaks  in  the  i'  echo  ani  $ is  an 

r Yr 

estimate  of  target  aspect  angle  on  the  r^  echo. 

Because  sub-targets  can  have  certain  range  rates  and  rate  of  change  of 
aspect  angle,  two  other  reasonable  features  are 


R 


\&t  / ( v— i ) L t L,- 1 

^ h JL 


• V— 1 * 

_i  y lirfilir 

\At/  (v-1)  L t t 
i*l  ^ r 


r 


(9) 


where  t ♦ is  the  time  between  pings. 

r*i  r 

A * 

Features  l,  a , N,  Q- f'j,  and  tf~)  are  significant  features  for  detecting 
submarine  targets.  To  achieve  a spherical  density,  v;e  vd.ll  need  the  following 
estimates  of  variarcss  of  these  features  for  sufc  targets: 

A A 

s£  - variance  (f)  for  sub-targets 

A * 

s ~ variance  (o)  for  sub-targets 

s 

s^j  ~ variance  (K)  for  sub-targets 

s .n  - variance  (4r)  for  sub-targets 

<£> 

A 

* A • 

s ..  *=  variance  (~-)  for  sub-targets 

<£>  6 

B.  Pesi^ni nr  Ron-linear  Filters  for  Potoctln/*.  JJonsnb_s_ 

Because  the  number  of  highlights  is  not  any  particular  fixed  value  for  many 
types  of  nonsubs,  a reasonable  additional  feature  for  nonsubs  is  the  average  change 
in  number  of  highlights, 

v— 1 

<»»>  * (^I)  l <Vl-  V (10) 

r~l 

Some  nonsubs  produce  multiple  targets  in  range.  Therefore  a reasonable 
feature  is  the  number  of  detected  targets  in  range,  TR,  an  estimate  for  which  is. 


v 

(TR)  *-  J Y (TU)j 


(ID 


r=l 


th 


vdt.h  (TR)^  the  number  of  targets  in  range  on  the  r echo  return. 

In  like  manner,  the  number  of  targets  in  bearing  TO,  is  a reasonable 
feature  for  certain  kinds  of  nonsubs  where  a reasonable  estimate  is 


(TO ) * y X (T0)i 
r=l 


(12) 


with  (TH)  the  number  of  targets  in  beaiirg  on  the  r ‘ echo  return. 

r 

For  certain  non  subs,  the  echo  length  my  change  such  that  a reasonable 
feature  is 


A 1 V'  ^ * 

Ae  a 7-T  1 Ivr  'r I 


where  is  the  e stirrated  echo  length  for  the  r echo. 


The  features  specifically  selected  for  nonsubs  thus  are  (a’>),  (TK),('19), 


and  (&o) . 
C.  Co'bir 


t_  .of.  Features 


The  vector  of  combined  features  for  subs  and  nonsubs  is  thus 

A A 

y = [i,  Or  N,  ) , (fy>'  (A‘0,  (Tit),  (TO ) , Ae]  (-U) 

where  the  ones  in  the  first  segment  of  the  partition  are  the  sub-features 
(denoted  v ) and  those  in  the  second  segment  are  the  nonsub  features  (denoted 

y ).  This  new  feature  vector  y has  features  (in  y ) specially  designed  for 

rs  a 

detecting  non-submarine  targets.  Because  there  undoubtly  are  correlations  between 
v and  v , it  is  not  wise  to  use  v as  the  feature  vector  unless  a priori 
knowledge  about  these  correlations  is  introduced.  Put  differently, ^ would 
contribute  nuiscance  features  when  trying  to  detect  a submarine.  In  the  following 
section  it  is  shown  how  a very  elementary  form  of  correlation  can  be  inserted  by 
using  conditional  density  functions. 

D«  Second  level  of  Correlation 

The  features  in  v (15)  may  not  have  spherical  statistics.  First  of  all, 

A A A A 

I and  a are  correlated  for  nor ubs . We  know  that  for  nonsubs,  both  t and  a 
can  be  larger  than  values  for  subs  or  both  can  be  smaller  than  valuer  for  subs. 
This  suggests  using  the  new  feature 


f [if  + (of 


for  nor.subs.  This  feature  might  also  be  satisfactory  for  subs.  Too  large  or 
too  small  values  would  signify  a sub. 


the  construction  cf 


certain  new  features  as  follows 


then  observe  the  chance  in  (c)  among  the  several 
targets.  A significant  change  indicates  a nonsubmarine 


target 


then  observe  the  pair  (X,  o)  pair  for  change.  Significant 


difference  of  this  pair  from  sub-values  indicates  a non 


submarine  target.  Also,  the  average  value  of  X over 


these  targets  can  be  used 


Another  nevx  feature  for  nonsubs  is 


oxpect  (a h)  to  be  increased  if  the  target  length  estimator 


because,  for  nonsub: 


X increase 


Introducing  Correl at ion  Lisina 


Functions 


Let  s index  the  sub-class  and  ns  index  the  nonsub  clas 


conditional  densities  of  can  be  written, 


Let  d(y  ) — i mean  that  the  decision  rale  using  feature  vector  y decides 

/*'S  ^-s 

class  i and  similarly  for  ^(yns).  Then  the  Bayes  framework  for  dimensionality 
reduction  suggests  the  following  approximations : 


_/  | \ _ /d(v  |nr.,d(y  ) — ns),  d(v  ) - 

P(y  ns)  - ■<  'Mas  1 * t-us  '*  V'di, 

P(^nSlns»d(^s)  = s)  > d(A,  > = 


ns 

s 


(16b) 


We  see  from  (16a)  how  the  non-sub  features  y contribute  in  a simple  v;ay  to 

#v  ns 

detecting  a submarine.  Eq. ’s  (16)  suggest  the  Block  diagram  shown  in  Figure  1 
where  all  the  nonlinear  a prion  knowledge  has  been  used  to  construct  the 
decision  rule. 

Because  nonlinear  functions  have  been  used  to  design  yg  and  yng  we  suspect 
that  the  statistics  of  the  class  conditional  densities  on  the  right  hands id e (16) 
are  spherical  and  unimodal.  Thus  v;e  represent  these  densities  with  Gaussian 
functions.  Define  the  mean  vectors. 


A>ss  = EUSIS>  ^ 

«®ss,ns  = EWS>  d(^nS)  = ns^ 

n *=  E[v  Ins,  d(y  ) - ns] 

*ais,ns,ns  L^-rs 1 J 

n = E[,v  Ins,  d(y  ) = ns] 

~ns,ns,s  v^ns 1 >^s  ’ 

Likewise,  define  the  corresponding  diagonal  covariance  matrices  T , T , 

✓>'Css  —S3, ns 


a?) 


I)  , and  v . Then, for  example, 

*hs,ns,ns  ~ns,r.s,s  * 

p(y  Is,  d(v  ) = s)  = expl"-g(v  - m )2fE  (y  - m )1  (IS) 

-^5  1 * >^S  jj.  L -*S  ~SOS  V-SSS  ~0  ^sss  J 

r\  A 1. 

where  5 is  the  dimensionality  of  y and  estimated  mean  vectors  and  cover-’ ance 
matrices  are  indicated. 

This  new  approach  to  pattern  recognition  is  quite  viable  because  we  are 

learning  more  about  sonar  every  day.  We  propose  to  continue  its  development  with 
the  objective  of  improving  performance. 


I -iuA 
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III.  Data  Ease 

The  data  ' ase  currently  being  used  is  fro::".  the  Rogers  Data  Base  obtained  iron 
ARL,  Austin,  Texas,  6/21/1SV0.  This  data  is  recorded  in  Figure  2.  Hie  PME  Reel 
inclusion  corresponds  to  the  corresponding  tape  at  Purdue.  Sore  of  the  measurements, 
speed,  class,  aspect,  and  number  of  pings,  are  included  in  this  table. 

After  using  this  Rogers  data,  we  expect  to  process  Sarsfeld  data  where  ranges 
are  between  2,000  yards  and  10,000  yards  with  some  as  far  as  40,000  yards  (subs 
only).  For  nonsubs,  the  maximum  range  is  7,000  yards.  All  the  Sarsfeld  data  is 
2 n.s.,  RDT,  c.w;  tracks  are  mostly  straight  line  but  some  are  circling. 

Another  source  of  data  is  Uitek  data  involving  ranges  out  of  20,000  yards. 

This  data  involves  very  large  sequencics  for  the  same  target  (200-300  pings),  for 
both  submarines  and  nonsubmarincs . The  Uitek  data  consists  of  close  to  45  reels 
of  good  data  - more  than  for  Sarsfeld. 

IV.  Performance  Using  Rogers  Data 

We  have  been  working  with  the  Regers  date  base  extracting  features  from  subs 

and  features  for  nonsubs.  Much  of  the  work  involves  finding  peaks  and  flat  places 

in  the  sonar  return.  Also,  the  target's  travel  is  fitted  with  a polynomial  to  get 

an  estimate  of  track  aspect.  Features  extracted  for  subs  are  UP  *»  number  of  peaks 

3/2 

NF  = number  of  flat  places,  = echo  effective  length,  ER  «=  (echo  energy)  (Range) 
3/2,  and  minimum  distance  between  the  largest  peak  and  its  nearest  peak.  A com- 

*J- 

puter  print  out  of  these  features  for  sub  echoes  and  nonsub  echoes  are  shown  in 
Figures  3a  and  3b.  In  Figure  3c,  track  aspect  angles  for  the  sub  targets  arc  shown. 
In  figure  4n  and  Figure  4b  is  shown  the  actual  computer  output  display  of  the  first 
two  sub  echoes  with  peaks  on  the  echo  marked  and  with  the  above  mentioned  features 
displayed.  In  Figures  4c  and  4d  we  see  computer  output  displays  for  the  first  two 
nonsub  echoes.  Figure  4e  shows  the  tracking  ships  track  and  the  target  track,  while 
Figure  4f  is  a blow-up  of  the  target  tra.ck. 


+ 


Sub  echoes  art  denoted  10F  and  noncub  echoes  1011  in  these  figures 


V.  Reverberation  Correlation  Cal c» lotion 3. 


■introduction 


. . . . pj 

Given  N sequential  waveforms  [v^(t ) j=i  > ^ro:n  M corresponding  band- 
pass stochastic  processes  obtained  from  M respective  channels,  it  is 

desired  to  compute  the  M2  possible  sample- ensemble-cross-co variance  functions 

{^mil(t'l>t2^n=i  d^i  defined  as  follows: 


rmn(tl’t2)  ‘ N 1 " vm(tl)][vsn(t2)  " %(t2)] 


where  v (t.  ) denotes  the  sample  mean  for  m 1 channel: 
m k 


%(tk}  ~ N IVsm(tk}  5 m = 


In  a digital  approach  a discretized  representation  of  the  waveform 
v^j(t^),  k = 1,2, ...,L  will  be  available.  If  the  stochastic  sampling  theorem 
is  utilized,  then  the  function  K^Ct-j^t^)  obtainable  from  signals  v^(t) 
sampled  at  a frequence  N less  than  the  My quest  frequency,  f = 2(fQ  + ^/2)j  fc 
is  the  sampling  frequency,  f is  the  bandpass  center  frequency,  and  W is  the 

low-pass  bandwidth  of  the  stochastic  processes  x^>  (containing  waveforms 

1 2 

y^j(t)  and  X^j(b)  respectively)  in  the  Rice  representation  of  "V  , i.c., 

vgj(t)  - xgV(t)  cos  iuot  - x<f(t)  sin  uuQt 

v €V  x(l)€v(l)  x(2)€v(2) 

*!T  vj’  sj  * » sj  * *j 


If  f » V/,  it  nay  be  practically  impossible  to  sample  at  the  resulting 
frequency  fg.  Moreover,  one  is  often  interested  in  computing  only  the  envelope 
of  ^(t^,  t^)  in  which  ca<5e  sampling  at  f may  be  inefficient.  An  alternate 
method  for  computing  this  envelope  denoted  ftCS^ftpt  ))  or  more  briefly 

Cmn^tl»t2^* 

There  are  M(m+l)/2  distinct  cro^s-correlation  functions. 


is  to  quadrature  sample  the  low  pass  functions  y.j^(t),  y ^ j (t > at  their 
Nyquist  frequency  of  W/2  hz.  One  may  then  compute  the  U?  low-pass  covariance 
functions  o^^i’  tj): 

N 

0k!.(ti’  \)>  * N IKmV  • ’ftV]  Lxm(tJ)  " ,n'(t3)] 
r=l 

k “ 1,  2;  l “ 1,  2;  n •=  1,  2,  M;  m «=  1,  2,  M (/*) 

It  is  well-known  that  for  jointly  stationary  stochastic  processes  , 

V , that  t (t,,  t„)  may  be  determined  from  the  /,  low-pass  functions 
m mn  1 < 

o-(tr  t^);  k *=  1,  2,;  l - 1,  2.  Specifically: 


°k!(tl’  t2)  “ °kV(V  tl»  0)  « “ 1.  2;  K - 1,  2 

°n^t2”  tl’  = °22^t2~  tl»  0^>  cr12^t2*’  tl»  ^ “ °2l(t2~  tl» 


emn(tl'  l2>  •’  i>n<V  tl-0)?+  L°"(V  ‘l-  0)? 


(5) 


Unfortunately,  Joint  stntlonarlty  Is  a condition  seldom  met  In  many  problems, 
e.g.  Sonar. 

B.  Procedure 

It  is  desired  to  sample  a bandpass  stochastic  process  at  the  lowest  possible 
frequency  which  allows  reconstruction,  one  such  technique  is  Quadrature  sampling. 
We  may  represent  the  sample  bandpass  function  x(t),  bandwith  W as: 

x(t)  *=  x^(t)  cos  WQt  - xQ(t)  sin  (i'ot 

I 


vhero  x^(t.)  and  xo(t)  arc  lowpass  stochastic  proccssos  of  which  the  highest 
frequency  is  \l/2.  Thus  we  my  apply  the  stochastic  sampling  theorem  to  the 
lowpass  proccssos  and  need  sample  x^(t)  and  xo(t)  at  W hz.  each.  Tho  envelope 
of  x(t)  is  Riven  l>y: 


e(x(t))  « Jk(( t)  -t  x*(t>  ' 

,pX  (t)  -1 

Tho  phaso  is  Riven  by  /(x(t))  ® fan  By  choooitiR  the  campling  frequency 

fg  to  bo  a submultiple  of  3 ffl  > W,  wo  may  obtain  the  samples  of  x^(t ) by 

simply  sampling  x(t)  at  frequency  f . by  taking  another  sample,  exactly  90° 

after  wo  take  each  sample  of  x.  (t),  wo  obtain  the  samples  for  xo(t).  Thun  wo 

seo  that  it  is  necessary  to  know  m and  to  generate  one  set  of  samples  at  f 
a 0 U'  6 

**  r~fr  hz  and  another  set  delayed  by  7^  radius.  Given  a synchonization  signal 

/* 

cos  o>  t,  it  is  desired  to  generate  a set.  of  sampling  pul  sen  at  a frequency 

°«<  O' 

f *=  N is  an  integer  and  another  set  delayed  by  t,  ! ~:Q  radius.  In  the 

8 U /» 

case  of  intorcst  to  *=  18  KHz,  N 18.  Then  f 1 1 Hz  and  we  may  adequately 
sample  a signal  with  a 1 KHz  bandwidth.  The  delay  between  signals  will  be 

td  A X 18  m:i0C*  72  nisec* 


In  tho  original  process,  these  samples  were 
sampling  at  /»  KHz  and  discarding  the  and  /;*h 


generated  erroneously  by 
samples  of  every  four  thus 


obtained. 


C.  Dcscr ipt Ion  of  Reverb  Correlation  Process 

Four  Gtuvoo  1,  AA,  A7,  AR  wore  quadrature  sampled  at  1 KJ!z,  starting  at 
2 sec.  after  tho  zero  time  pulse  to  obtain  Cl^  in  phase  and  6/*  quadrature 
samples  for  each  stave.  A total  of  1 80  such  sets  of  samples  (A  x 180  *-  7.70 

■t 

in  phase  vectors  and  720  quadrature  vectors)  we  obtain  and  stored  an  IBM 
1130  Magnetic  Disk. 

The  correlation  coefficient  functions  we  computed  for  four  different 
subsets  of  tho  180  vectors  (#  1-60,  H 60-121,  #121-180,  and  H 1-180)  in  order 
to  examine  tho  nonstationarity  of  tho  process.  For  each  such  ensemble,  the* 

66  possible  distinct  correlation  coefficients  were  computed  as  follows.  For 
each  pair  of  the  four  staves,  there  are  two  in  phase  and  two  quadrature  vectors. 
For  each  of  these  four  vectors,  the  means  wore  computed  and  then  subtracted  out. 
The  four  6A  x 6A  covariance  matrices  wore  then  computed.  These  were  then  com- 
bined using  the  relation  derived  earlier  to  obtain  the  sampled  envelope  of 
the  cross  covariance  function  for  that  pair  of  staves  t! (!I1)(n),  tj.^0  COJ*reia- 


mn 


tion  coefficient  matrix  p was  then  computed  using  tho  relation: 


mn 

liJL 


nn  _ 

Vc,ij  °u 

D.  Experimental  Stave  Correlations 

We  were  asked  by  Naval  Ships  System  Command  to  investigate  correlations 
among  multistave  data.  The  procedure  for  numbering  staves  is  shown  in  figure  5. 

Echos  or  pings  were  from  3 different  modes.  The  ping  numbers  and  correspond- 
ing modes  are  shown  below: 

Ping  # Pulse  Width  (at  71;  ips) 

1 - A8  60  m.  sec.  e.w. 

49-100  Chatterbox  (^pulses,  each~’10m.  sec.) 

101-180  Mixed  e.w.  (A  short  followed  by  two  60  m.scc.) 


Signals  for  the  various  staves  were  recorded  on  14  channel  tape  as  illustrated 
In  Figure  5.  A general  block  diagram  of  the  procedure  Is  shown  in  Figure  7. 

A sequence  of  pictures  shown  In  Figure  8 show  the  reverberation  correlations  c- 
cross  staves  for  several  configuration.  A sequence  of  such  pictures  are  In  the 
procession  of  Navel  Ships  System  Command. 


Part  II 


FORWARD  (Next  90  Days) 


We  are  processing  more  of  the  Rogers  data  so  as  to  report  more  per- 


formance results.  We  will  derive  features  for  the  nonsubmarine  target 


such  as  sea  weed  where  there  are  multiple  echoes  in  either  range,  bearing. 


or  both  range  and  bearing.  We  are  inserting  a priori  knowledge  concernin'; 


existence  of  two  major  peaks  for  subs  and  peak  signal  power/average  noise. 


Sometime  soon  ue  must  begin  to  give  attention  to  "false  targets"  and  video 


scan  data  to  achieve  a more  accurate  estimate  of  target  bearing.  This 


estimate  of  target  bearing  is  important  when  estimating  target  spaced 


and  target  length. 


i i 


Figure  2:  ROGERS  DATA  Austin  6/21/70 
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Figure  2 continued 
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Figure  4a:  Sub  Echo,  Features,  and  Marked  Peaks 
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Figure  4b:  Sub  Echo 
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Figure  4c:  Nonsub  Echo 


RELATIVE  STAVE  LOCATIONS 
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Figure  5:  STAVE  LAYOUT  FOR  TAPE  # RVB-1 
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Direct 

Stave  //I 
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Direct 

Stave  ft 2 

3 

Direct 

Stave  //4” 

4 

Direct 

Stave  #3 

5 

Direct 

Stave  #47 
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Direct 

Stave  #4 
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Direct 

Stave  #46 

8 

Direct 

Stave  #5 

9 

Direct 

Stave  #45 

10 

Direct 

Stave  #6 

11 

Direct 

Stave  i! 44 

12 

FM 

12.95  KHZ  reference 

13 

Direct 

Stave  #43 

14 

FM 

1 sec  pulse;  starts 

when  bean  is  dead  ahead 

Recorded  C:  7^  ips. 

Set  footage  counter  to  zero  0 first  trace  of  audio  on  stave  channels; 
then 

ft-,  Pin  ns 

0-438  1-48  Single  ping  CVJ,  Long  pulse  (60  m.scc) 

456-938  49-100  Chatterbox  ( 30  closely  spaced  pings) 

953-1505  101-178  13  burst  of  4 long,  two  short  CW 


Figure  6 : 


MULTI-STAVE  REVERR'C RATION  DATA 


